Fractional Klein-Kramers equation for superdiffusive transport: normal versus 
anomalous time evolution in a differential Levy walk model. 
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We introduce a fractional Klein-Kramers equation which describes sub-ballistic superdiffusion in 
phase space in the presence of a space-dependent external force field. This equation defines the 
differential Levy walk model whose solution is shown to be non-negative. In the velocity coordinate, 
the probability density relaxes in Mittag-Leffler fashion towards the Maxwell distribution whereas in 
the space coordinate, no stationary solution exists and the temporal evolution of moments exhibits 
a competition between Brownian and anomalous contributions. 
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Classically, Brownian stochastic transport processes 
in the phase space spanned by velocity v and coordi- 
nate X are described by the deterministic Klein-Kramers 
equation (KKE). In the low and high friction limits, the 
KKE reduces to the Rayleigh equation which describes 
the relaxation of the velocity probability density function 
(pdf) towards the Maxwell distribution, and the Fokker- 
Planck-Smoluchowski equation controlling the temporal 
approach of the Gibbs-Boltzmann equilibrium, respec- 
tively H, E|. The KKE is therefore a fundamental equa- 
tion in non-equilibrium systems dynamics. 

Brownian transport is characterised through the Gaus- 
sian pdf and the linear time dependence (a;^(t)) = 2Kt 
of the mean squared displacement in the force-free diffu- 
sion limit, its universality being guaranteed by the central 
limit theorem 0. In a broad variety of systems, however, 
it has been found that correlations in space or time give 
birth to anomalous transport whose pdf is non-Gaussian 
and/or whose mean squared displacement is non-linear 
in time ||3|. These systems include charge carrier trans- 
port in amorphous semiconductors ||4], tracer dispersion 
in convection rolls and rotating flows ||5|, |6[ , capillary sur- 
face waves [0, the motion of bacterias and the flight of 
an albatross S, intracellular transport [B, transport in 
micelles [[l^l, 2-D dusty plasmas [^, the dynamics in 
(bio) polymeric systems |1^, and the NMR diffusometry 
in porous glasses and percolation clusters [|l3[, among 



others. 

The systems we are interested in fall into the broad 
class whose force-free diffusion behaviour is characterised 
through the power-law form (x^(i)) oc t^ , which sepa- 
rates into subdiffusion (0 < x < 1) and superdiffusion 
{k > 1). The continuous time random walk (CTRW) 
model has proved to be a well-suited framework which 
accounts for such anomalous diffusion for the entire spec- 
trum of X jlj]. Especially in the sub-ballistic superdif- 
fusive domain 1 < >f < 2, Levy walks which couple long 
ffight times with a time cost have been a successful tool, 
e.g., in fluid dynamics Iq, na. The space-time coupling of 
Levy walks leads to finite moments and they are therefore 
fundamentally different from Levy flights which exhibit 
a diverging variance pq | . 

In the presence of external force fields, the CTRW ap- 
proach is less flexible. It has been realised that fractional 
equations constitute a tailor-made framework to formu- 
late the underlying dynamics equations in coordinate and 
phase space; see, for instance, |1^, |l^, |8[ |l^, g^l ^^^^ ref- 
erences therein. In the subdiffusive domain < >f < 1, 
a fractional KKE was derived from CTRW models |2l|] 
and from the Chapman-Kolmogorov equation |g^. A 
similar consistent generalisation to systems in the regime 
1 < X < 2 is still outstanding. In this note, we propose 
the fractional KKE 
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for the description of sub-ballistic superdiffusive anoma- 
lous transport. In Eq. (|l|), P{x,v,t)dxdv is the joint 



probability to find the test particle of mass m with co- 
ordinate X, . . . ,x + dx and velocity v, . . . ,v + dv at time 



t. 7 denotes the friction constant which quantifies the 
effective dissipative interaction with the environment, k 
is the velocity diffusion constant and F{x) — —d^{x)/dx 
is an external force field. The factor g^ has dimension 
\Qa\ = sec~" and is a function of some time scale t char- 
acteristic for the waiting time process and of an inter- 
action time scale t* [p2|, see also below. In the limit 
a = 1, Eq. (nl) corresponds to the standard KKE whereas 
for < a < 1, the process contains a scale- free mem- 
ory of the power- law type entering through the fractional 
Riemann-Liouville operator qD\~°' = d/dti^aD^"^ with 
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dt'W{v,t'){t-t')'^-^/T{a). (2) 



The operator o^T" possesses the important property 
/o e~"* o-CT"/!^) = w^"/(u) under Laplace transfor- 
mation. 

With respect to its coordinate x, Eq. (|l]) is a Liouville- 
type equation which gives rise to the relation ■^{{x{t))) — 
{{v{t))) between the velocity- and coordinate-averaged 
random variables x and v; by this, the KKE (mj differs 
from the subdiffusive model in Ref. [E2|. The action of 
the friction, force and velocity diffusion in Eq. p) enters 
through the non-local memory relation brought about by 
the fractional operator; i.e., the local fluctuations of ve- 
locity which follow from a Fokker-Planck-like equation 
in the Brownian case are now described by a fractional 
Fokker-Planck equation. The connection to the space co- 
ordinate is assumed to be given through the classical drift 
term —dvP/dx. From these basic physical assumptions, 
the fractional KKE (|l|) emerges. 

Let us explore the behaviour predicted by Eq. (|l|) in 
more detail. By integration over the spatial coordinate 
X, the fractional equation for the velocity pdf P{v, t) 
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obtains which reduces to its Brownian counterpart in the 
limit a = 1. 

In the force- free limit, (||) corresponds to the fractional 
Rayleigh equation discussed in Refs. |2|, |2J]. The relax- 
ation of velocity moments can be inferred directly from 
Eq. (^ in the force-free case. Accordingly, the fractional 
relaxation equation |Tq, Oa] 

i{v{t)) = -^o,oDl~''{v{t)) (4) 

emerges where ja = Qal- For vq — w(0), its solution 



(w(i)) = voE^ {--tct") 



(5) 



features the Mittag-Leffler function Ea (— 7ai") = 
S^o (~7"^")" /^(-'^ + Q^") which is monotonically 
decaying and which interpolates between the ini- 
tial stretched exponential behaviour Ea {—Jat") ~ 



exp (— 7Q,i"/r(l + a)) and the final inverse power-law 
pattern E^ {--fat°') ~ (7ai°T(l - a))~^ (26|. Thus, the 



mean relaxation time diverges. The second moment 

(v^t)) = Kh + {vl - Kh)Ea (-27„r) 



(6) 



equilibrates Mittag-Leffler fashion towards the equilib- 
rium value k/7, and the pdf P{v,t) relaxes towards a 
Gaussian. In thermodynamical systems, by compari- 
son, the velocity pdf follows the Maxwell distribution 
P{v) = (27rfcBr/TO)"^/^exp(-TOuV[2fcBT]V so that we 
find the Einstein relation n = kBTj/m [n6l O. The 
velocity-velocity correlation function associated with Eq. 
(|)iB 
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whose long time behaviour is equivalent to the CTRW- 
Levy walk result where the t~"-scaling follows from the 
behaviour of the cumulative waiting time distribution. 

This corresponds to the physical model underlying Eq. 
(|l|): in the velocity space, we assume that the particle 
undergoes collisions, i.e., random velocity changes, such 
that successive collisions are separated by time spans gov- 
erned through the long-tailed waiting time pdf V'(t) ~ 
Aa{i/T)^^^°' whose characteristic time T = J^ i'ij}{i)di 
diverges. In this case, the velocity-velocity correlation 
function is proportional to the probability that no scat- 
tering took place before time i, and thus {v{Q)v{t)) — 
Vq J^ ip{t)dt oc z)g(i/r)~". In addition, we assume that 
the external force field F{x) constantly acts on the test 
particle during its sojourns. 

The characteristic behaviour of this model is that it 
includes arbitrarily long sojourns but leads to finite mo- 
ments of any order and an exponentially decaying pdf 
P. We call the process associated with the fractional 
KKE (|^) a differential Levy walk model noting that it 
is an approximation to CTRW-Levy walks which offers 
the distinct possibility to study the effects of an external 
force field on the behaviour of lower order moments of a 
Levy walk. 

The velocity average of the fractional KKE (|^) can 
be performed by integration over J dv and J vdv, and 
combination of the two resulting equations. With (w^) — 
ksT/m, this procedure yields the fractional telegrapher's 
type equation 
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which exhibits a transition from a short-time ballistic 
behaviour with {x'^{t)) ^ {K^)t^ in the force-free case, 
to the long-time or high-friction limit governed through 
the fractional Fokker-Planck-Smoluchowski equation 
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Its force-free limit exactly corresponds to the enhanced 
diffusion approach reported in Ref. |£^ which exhibits 
a distinct bimodal behaviour characterised through two 
separating humps. From Eq. (S) one observes the differ- 
ence to the fractional Fokker-Planck equation for sub- 
diffusive processes from Ref. 
temporally local with dP/dt. 
K = ksT / {m'^a)', note that this does not represent a 
real Einstein relation as Eq. ffl) does not have a station- 
ary solution, see below. The mean squared displacement 
corresponding to Eq. (0) with F{x) = is given through 



20 : here, the force is 
In Eq. (0), we identify 
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which describes sub-ballistic superdiffusion. This result 
is equivalent to the fractional Kramers model discussed 
in Ref. ||, and to the CTRW-Levy walk Q. It is in 
the presence of an external force that the fractional KKE 
(|l|) exhibits a fundamentally different behaviour. 

To see this, consider the fractional Fokker-Planck equa- 
tion (H) for non-trivial types of the external force F{x). 
Accordingly, the first moment is given by -^{x{t)) — 
{F{x))/{m'y) which can be solved for constant or linear 
forces. In particular, for the constant drift F{x) — Vm^, 
the first moment becomes 



{x{t)) = Vt 



(11) 



which corresponds to the traditional drift behaviour pj. 
The second moment 
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combines this drift with the sub-ballistic behaviour ex 
t2~« such that the variance {Ax{t)'^) = {{x^ [t)) - {x{t))'^) 
is given by Eq. (p^). This behaviour is analogous to 
the Galilei invariant diffusion-advection model derived 
in Ref. pq] for the subdiffusive case. In our case, the 
analytical solution is given by the free (superdiffusive) 
solution explored in detail in Ref. [ETI, taken at the 
translated coordinate x — Vt. Here, the two-hump so- 
lution travels with velocity V and is symmetric to the 
point X{t) = Vt. It should be noted that due to the 
behaviour of the first and second moments, a connection 
of the form {x{t))v oc V{x^{t))a (see, e.g., Ref. 0, |l7|) 
between the first moment in the presence of the constant 
drift V to the second moment in absence of this drift, 
does not exist, in contrast to the subdiffusive case |2(](| . 

Similarly, for the Ornstein-Uhlenbeck potential $(x) — 
^muj'^x'^ which exerts the linear restoring force F{x) = 
muj'^x, the first moment shows the exponential relaxation 
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which contrasts the Mittag-Leffler patterns recovered in 
the subdiffusive model in Refs. p3| as well as the su- 
perdiffusive fractional KKE presented in Ref. p3]. The 



second moment has the Laplace transform 
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whose inversion leads to 
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which combines the exponential relaxation of the ini- 
tial condition, which was already found for the first 
moment, with the confluent hypergeometric function 
iFi. Note that the second term in Eq. (|l5|) is 
equal to the expression 2Kt^~°'Ei^3-a (— Sw^i/TJ = 
2^^-" ^~^o (-2^2^/7) /r(3 - a -I- n) in which we used 
the generalised Mittag-Leffler function £'i.3_q, |2g]. The 
variance is consequently given in terms of 



Ax{t)^ = 2Kt^-"Ei^, 
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which interpolates between the freely diffusive behaviour 
(|1Q| ) for short times, and the long time power-law pattern 
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Accordingly, the mean squared displacement Ax{t)'^ in- 
creases in the course of time even for long times. It is 
straightforward to show in general that the solution of 
the fractional Fokker-Planck-Smoluchowski equation (0) 
does not have a stationary solution. This is exactly what 
is expected from a Levy walk whose main characteristic 
is the continuous approximation to a Levy stable distri- 
bution |I|]. 

Let us now address the problem of negativity p9| which 
was encountered for the fractional KKE proposed in Ref. 
Q in the overdamped limit. Firstly, it is possible to 
show that for all constant or linear forces, moments of 
any order are positive. This is due to the fact that our 
model equation (|l|) leads to (generalised) Mittag-Leffler- 
type relaxation patterns which are known to be strictly 
monotonically decreasing and which are everywhere pos- 
itive Ig^. In contrast, the model from Ref. ||2J] leads to 
Mittag-Leffler indices between 1 and 2 and therefore to 
temporal oscillations which in turn give rise to negative 
portions in the "pdf" for shorter times [g9|. The posi- 
tivity of the solution of our fractional KKE (||) can be 
demonstrated also for an arbitrary external force fleld. 
We know that the velocity pdf P{v,t) is positive for all 
times [^ Q, and that the position-space pdf P{x,t) is 
positive in the force- free case [E^. As the force enters 
through the local drift term, also the pdf in the pres- 
ence of an arbitrary external force F(x) is positive, and 
consequently the joint pdf P{x, v, t) is positive. 

Formally, the fractional KKE (Q) corresponds to the 
"Levy rambling" model developed in Ref. [p^ but with 
the force entering symmetrically to friction and diffusion 



(i.e., {v) = -r]VT* + F{x)t* /m in Eq. (62) of Ref. H). 
Thus, Eq. (P corresponds to an approximation to the 
CTRW-Levy walk model in which the particle under- 
goes continuous motion between scattering events which 
change the velocity coordinate. Friction and force re- 
main events which enter through an effective time scale 
T*, i.e., they correspond to point-like interactions in the 
long time limit t ^ max{r, t*}. This gives rise to the 
fact that the spatial distribution does not reach an equi- 
librium state, or, in other words, that space and time do 
not decouple in the underlying equations (if) and (pf) so 
that the variables x and t cannot be separated, also a 
typical property of CTRW-Levy walks. 

We have introduced a new fractional approach to the 
phase space description of superdiffusive sub-ballistic 
transport processes. The obtained fractional KKE leads 
to the Mittag-LefHer relaxation of the velocity distribu- 
tion towards the classical Maxwell-Boltzmann equilib- 
rium. In contrast, the long-time or high- friction limit 
of the spatial distribution does not possess a stationary 
solution. The process is characterised by a combination 
of the classical time evolution found for the first mo- 
ment, such as the exponential relaxation of the initial 
condition in the presence of a linear force field, with a 
time-dependence which is governed by the fractional or- 
der a, i.e., the "memory strength". The fractional KKE 
fulfils a generalised Einstein relation in velocity space. 
In coordinate space, no analogous generalisation of the 
Einstein relation exists. These properties set the present 
approach apart from previous fractional models and we 
therefore regard it a distinguished candidate for the de- 
scription of enhanced, sub-ballistic transport in external 
force fields, and as complementary to the established sub- 
diffusive fractional KKE-description. 

We call the associated stochastic process a differential 
Levy walk in reference to the fact that Eq. (|l|) correctly 
describes lower order moments of a CTRW-Levy walk in 
an external force field and in phase space. A similar in- 
vestigation of such situations within the original CTRW 
scheme would be very tedious. We expect that this study 
will contribute to a better understanding of superdiffu- 
sive transport and its dynamical foundation. 
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